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Contamination on Lyman continuum emission at z ^ 3: 
implication on the ionising radiation evolution 
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ABSTRACT 

We investigate the possibility of contamination by lower-redshift interlopers in the measure 
of the ionising radiation escaping from high redshift galaxies. Taking advantage of the new 
ultra-deep VLT/VIMOS U-band number counts in the GOODS-South field, we calculate the 
expected probability of contamination by low redshift interlopers as a function of the U mag- 
nitude and the image spatial resolution (PSF). Assuming that ground-based imaging or spec- 
troscopy can not resolve objects lying within a 0.5" radius of each other, then each z>3 galaxy 
has a 2.1 and 3.2% chance of foreground contamination, adopting surface density U-band 
number counts down to 27.5 and 28.5, respectively. Those probabilities increase to 8.5 and 
12.6%, assuming 1.0" radius. If applied to the estimates reported in the literature at redshift ~ 
3 for which a Lyman continuum has been observed directly, the probability that at least 1/3 of 
them are affected by foreground contamination is larger than 50%. From a Monte-Carlo sim- 
ulation we estimate the median integrated contribution of foreground sources to the Lyman 
continuum flux (/900). Current estimations from stacked data are >2<x of the median inte- 
grated pollution by foreground sources. If the correction to the observed /900 flux is applied, 
the relative escape fraction decreases by a factor of ~ 1.3 and 2, depending on the cases re- 
ported in literature. The spatial cross-correlation between the U-band ultra-deep catalog and a 
sample of galaxies at z>3.4 in the GOODS-South field, produces a number of U-band detected 
systems fully consistent with the expected superposition statistics. Indeed, each of them shows 
the presence of at least one offset contaminant in the ACS images. An exemplary case of a 
foreground contamination in the Hubble Ultra Deep Field at redshift 3.797 by a foreground 
blue compact source (U=28.63±0.2) is reported; if observed with a low resolution image (see- 
ing larger than 0.5") the polluting source would mimic an observed (f 1500/ f 900) obs ~ 38, 
erroneously ascribed to the source at higher redshift. 

Key words: Galaxies: evolution - Galaxies: high-redshift - Individual source: 
GDS J033236.83-274558.0 



1 INTRODUCTION 

Ionising radiation from star-forming galaxies is a plausible primary 
source of cosmic reionisation. The amount, the evolution with red- 
shift and its contribution to the cosmic reionisation is still a mat- 
ter of debate (e.g., Steidel et al. (2001), Giallongo et al. (2002), 
Shapley et al. (2006) [S06, hereafter], Siana et al. (2007, 2010), 
Faucher-Giguere et al. (2008), Iwata et al. (2009) [109, hereafter]). 

Malkan et al. (2003) and Siana et al. (2007, 2010) have stacked 
tens of deep ultraviolet imaged galaxies at redshift ~ 1 and no de- 
tection has been reported. Similarly, Cowie et al. (2009) have com- 
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bined ~ 600 galaxies observed with GALEX at the same redshift 
and also in this case the result was a non detection. At higher red- 
shift (z>3) an escape of Lyman continuum photons (LyC hereafter) 
seems to be detected in ~ 10% of the galaxies (e.g. 109 and S06), 
suggesting a possible evolution of the escape fraction with redshift. 
It is still unclear which is the typical escape fraction to be expected 
from starburst galaxies at redshift z>3. While some results of nu- 
merical simulations suggest low values (less than 10%, e.g. Gnedin 
et al. 2008), others claim escape fraction as high as 80% (e.g., Ya- 
jima et al. 2009, Razoumov & Sommer-Larsen 2009). The two pos- 
sibilities imply two completely different scenarios: in the first case 
high-redshift galaxies turn out to be inefficient in releasing ionising 
radiation into the intergalactic medium, conversely in the second 
case, they may play a dominant role in this redshift regime. 

As redshift increases, the probability of foreground contami- 
nation by blue galaxies that mimic the ionising radiation also in- 
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Table 1. Number counts derived from the ultradeep VLT/VIMOS catalog 
down to U ~ 30 (lcr limit within 1" aperture radius, Nonino et al. 2009). 
Completeness correction has been roughly applied at magnitudes fainter 
than 28.5 (~3cr limit). The values in parenthesis have been calculated by ex- 
trapolating linearly the number counts up to U-mag 30.5 by fitting a straight 
line in the range U-mag 25-28 with a slope 0.2 in the [LogN] vs. [U-mag] 
plan. 



U-ma^ 




Ncounts 


Cumul 


Ulow 


-Uup 


[deg~ 2 ] 


(m < Uup) [deg~ 2 ] 


22.0- 


22.5 


1400±40 


1400 


22.5- 


23.0 


2200+50 


3600 


23.0- 


23.5 


4400±70 


8000 


23.5- 


24.0 


8300±90 


16300 


24.0- 


24.5 


14900±120 


31200 


24.5- 


25.0 


23800±150 


55000 


25.0- 


25.5 


33500±180 


88500 


25.5- 


26.0 


45600±210 


134000 


26.0- 


26.5 


57900±240 


191900 


26.5- 


27.0 


72300±270 


264200 


27.0- 


27.5 


85300±270 


349500 


27.5- 


28.0 


90200+300 


439700 


28.0- 


28.5 


80400+280 


520200 


28.5- 


29.0 


54900+230 (177800) 


575100 (698000) 


29.0- 


29.5 


23700+150 (223900) 


598800 (921900) 


29.5- 


30.0 


4500+70 (281800) 


603300 (1203700) 


30.0- 


30.5 


— (354800) 


603300 (1558500) 



creases. Therefore it is important to study the contaminated sys- 
tems and try to quantify their occurrence. The recent ultra deep 
VLT/VIMOS U-band observations in the GOODS-S field (Nonino 
et al. 2009) allow to investigate this issue to unprecedented depth 
and to better calculate the probability of such contamination (a 
comparison between the observed superpositions in the GOODS - 
S area and those expected is also performed). Very deep expo- 
sures with exquisite spatial resolution, as the Hubble Ultra Deep 
Field (HUDF, Beckwith et al. 2006), make it possible to verify 
the results of such a calculation over limited areas of the sky. The 
work is structured as follows: in Sect. 2.1 the probability distribu- 
tions of foreground superposition are calculated and in Sect. 2.2 
it is compared to the observed direct detections reported in liter- 
ature. In Sect. 2.3 the Monte-Carlo simulations are described to 
derive the integrated contribution of foreground sources. In Sect. 
3 the U-band catalog is cross-matched to the redshift ~ 4 sam- 
ple and checked for contamination. Sect. 3.1 discuss the case of 
an exemplary false LyC escape detection in the HUDF. Sect. 4 
summarise the results. The standard cosmology is adopted {H -10 
km s~ l Mpc~ l , £2^=0.3, £2 A =0.7). If not specified, magnitudes are 
given in AB system. 



2 LYMAN CONTINUUM ESCAPE RADIATION AT Z>3 

2.1 The expected probability of foreground contamination 

Taking advantage of the ultra-deep VIMOS U-band imaging in 
the GOODS-S field (Nonino et al. 2009), that reaches magnitude 
29.8 (28.5) at lcr (3cr) within 1" aperture (radius), we determine 
the likelihood of foreground contamination of z>3, as these galax- 
ies must reside in the foreground and are emitting at wavelengths 
that mimic LyC of LBGs. Similar to the discussion in Siana et al. 
(2007) we assume that the spatial distribution of the foreground 
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Figure 1. Probabilities (Prob, Y-axis) that a fraction (%, X-axis) of an ob- 
served sample is contaminated by foreground sources. Calculations have 
been performed varying the single case probability, p, in turns related to 
the U-band number counts (from left to right panels) and seeing (differ- 
ent curves in each panel). Top panels: The cumulative probability distribu- 
tion versus fractions (%) of the entire sample (N). It is the probability that 
at least a fraction K(%) of the Af galaxies observed (P(^K) with N > K) 
is contaminated by foreground sources. In the fainter limit case, U-mag < 
29.5 (left panel), calculations have been performed with the observed num- 
ber counts (solid lines) and those corrected by uncompleteness in the range 
28.5<U-mag<29.5 (dotted lines, see text for details). The arrows show how 
the curves change from uncorrected to corrected. Bottom panels: The prob- 
ability distributions that a fraction K(%) of the N galaxies observed (P(K)) 
are contaminated. Increasing the U-band surface density and the seeing, the 
probability of superposition increases. In all panels the two vertical dotted 
(red) and dashed (blue) lines mark the percentages of the direct LyC detec- 
tions by 109 (17/198) and S06 (2/14), respectively. 

(U-detected) galaxies is uniform and uncorrelated with those at 
higher redshift. If we adopt the surface density of objects with 
U(AB)<28.5 (520200 deg' 2 , see Tableland Nonino et al. (2009)) 
and assume that ground-based imaging or spectroscopy can not re- 
solve objects which lie within a 0.5 (0.8)" radius of each other, than 
we would expect that each z>3 galaxy has a /?~3.2 (8.1)% chance 
of foreground contamination. Figure [T] summarise the probabilities 
that a certain fraction (%) of the total sample N is contaminated by 
foreground sources, as a function of the seeing and U-band num- 
ber counts. Q It is evident that increasing the U-band magnitude 
limit of galaxy counts and the seeing value, the probability of fore- 
ground contamination increases. Depending on the seeing condi- 
tions and the U-band magnitude limit adopted, the fraction of the 
population contaminated by foregorund sources ranges between 3 
and 15%, more severely if the faint fluxes are investigated (see Fig- 
ure[T|). Currently, the direct LyC detections reported in literature re- 
gard fractions smaller than 15%. It is worth noting that stacking of 
sources can also be affected by contaminantion, being a result of a 

1 The probability to observe K (or > K) contaminated sources in a sample 
of N (> K) high-z galaxies given the probability p of the single case is: 

f(K) = (%)p k (1- P f- K ; P(> K) = zl k /(0; 
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sum that may include contaminated cases not detected individually 
(see Sect. 2.3). In the following sections we consider the observed 
number counts down to U-mag 28.5 (3<x limit), beyond this limit a 
linear extrapolation is performed (see Sect. 2.3). In the next section 
the derived probabilities are compared to the observations of direct 
LyC. 

2.2 The current LyC detections 

S06 reported on LyC detections: 2 out of 14 galaxies show a clear 
signal below 912A rest-frame in their deep spectra. In particular 
SSA22a-C49 fe=3.115) and SSA22a-D3 (z=3.067) show a LyC 
flux of / 9 oo=(lL8+l.l)xlO- 31 ergs s~ l cm' 2 Hz~ l (mag(AB) = 
26.22) and / 9 oo=(6.9±1.0)xlO- 31 ergs s~ l cm' 2 Hz' 1 (mag(AB) 
= 26.80), respectively. Adopting the seeing value reported in S06, 
never worse than 1.0" and larger than 0.8", and the number counts 
in the magnitude range 26<U<27 (130200 deg~ 2 ), the probabilities 
that at least one galaxy (out of 14) in the S06 study is subject to 
foreground contamination assuming seeing 0.8, 0.9 and 1.0", are 
25, 30, 36% respectively. There is a 3, 5, 7% chance that at least 
two detections are contaminated. 

It is worth noting that in the case of spectroscopic observa- 
tions, the slit width (1.2" in S06) and its orientation on the sky may 
influence the net contribution of a close ("sub-seeing") contami- 
nant, since it may be partially included in the slit if there is a spa- 
tial offset. Therefore the /900 flux due to offset foreground sources 
is decreased. In principle, there is also the additional possibility 
to identify the foreground source through its spectral fatures if the 
emission lines are present in the wavelength range and are intense 
enough to be detected. 

109 reported similar LyC detections in LBGs and LAEs at red- 
shift 3.1 with dedicated narrow band imaging (NB359). In partic- 
ular 17 (7 LBGs and 10 LAEs) sources out of 198 (73 LBGs and 
125 LAEs) show a clear signal (>3cr). Looking at their Figure 4 the 
NB359 magnitude of the 17 cases is in the range 26-27.5, the sur- 
face density of U-band detected galaxies in this interval is 215500 
deg~ 2 (Nonino et al. (2009)). Assuming as in 109 that a foreground 
object within 1.0" radius cannot be distinguished in the image, 
each object has 5.2% chance of contamination, corresponding to a 
59(46)% chance that at least (10/11) sources out of 198 are contam- 
inated by foregrounds. If we adopt a U magnitude range 26.5-27.5 
as in 109, the U surface density becomes 157600 deg~ 2 , correspond- 
ing to a 49% of probability that at least 8 objects out of 198 are con- 
taminated. Furthermore, 109 note that shapes and positions of the 
emitting regions in the sub- sample of LBGs is different from those 
in the R-band, reporting an average offset of 0.97". Since spatial 
offsets are detected, 109 used adaptive apertures (centered in the R 
image) to measure the bulk of the fluxes in the R and NB359 filters 
and derive the colour (i.e. the observed /1500//900 ratios). The 
aperture diameters adopted vary from 2.0" to 4.0", therefore, apart 
from the "intrinsic" confusion due to the seeing, the relatively large 
apertures may include even more contaminants. Considering the 
sub-sample of 73 LBGs for which an offset has been measured, and 
adopting the median case of 1.0" radius and surface number den- 
sity of galaxies with 26.5 < U < 27.5 of 157600 deg~ 2 , it turns out 
that there is 53 (2) % chance that at least 3 (7) estimated colours are 
contaminated. With 1.5" radius the chance is 96 (44)%. 109 suggest 
that the spatial offset may be related to the intrinsic nature and/or 
the geometry of the escape radiation process, mantaining however 
open the possibility of a superposition for part of the sample. 

Finally, it is worth noting that the two LyC detected galaxies 
by S06 (SSA22a-D3 and SSA22a-C49) have been observed also 
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Figure 2. Upper panels: Simulations on the LyC contamination in the case 
of S06. In the left side the flux distribution of the 10000 rectangular ape- 
tures randomly positioned over the observed U-band image (dotted line) 
and the noiseless image (with number counts down to -28.5, solid line) are 
shown. As expected the noiseless image is peaked at the null flux value (cor- 
responding to rectangular apertures not intercepting sources). In both cases 
the asymmetry towards brighter fluxes is the signature of pollution by inter- 
lopers (the case with Umax=26 is shown). In the rigth panel the distribu- 
tions of the averages over 14 rectangular apertures drawn from the noiseless 
image are shown, in the case of Umax=24, 25 and 26 (dotted, dashed and 
solid lines, respectively). The vertical dotted lines from left to right mark 
the median, the 1 and 2 cr confidence levels in the case Umax=26. The es- 
timation reported by S06 (average of 14 spectra) is also shown with a filled 
circle. Lower panels: The same as upper panels, but the simulations have 
been computed for the case of Steidel et al. (2001). 

by 109. In the case of SSA22a-D3, 109 notice that the object is not 
visible in their NB359 image, even though the flux density limit is 
well below (two times) the limit reported by S06. The other source, 
SSA22a-C49, shows a NB359 signal (2.95cr level) spatially shifted 
with respect to the R band. From the ACS/F814 imaging (shown in 
109, Figure 2, top right panel) it is apparent a distinct substructure 
within 1.0" from the brighter LBG, spatially consistent with the 
emission in the NB359 image. It would be important to obtain the 
redshift of such companion. 

All these considerations suggest that the situation is still un- 
clear and even though the pure foreground contamination may not 
entirely explain the direct LyC detections at high redshift, its effect 
is not negligible. 

2.3 The integrated contribution of foreground sources 

In the case of the stacking of LyC non-detections the foreground 
contamination may play a role in the result. The stacking process 
increases the signal to noise ratio of the possible true emission of 
LyC, but at the same time a bias by foreground faint sources (not 
detected singularly) may arise. We calculate the expected average 
integrated contribution of the foreground blue sources in different 
rectangular apertures performing Monte-Carlo simulations. In or- 
der to avoid possible contribution of the error fluctuation of the 
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background, a noiseless image has been generated. It has been 
obtained by running the SExtractor algorithm (Bertin & Arouts 
1996) allowing the detection down to U-mag^28.5 and requiring 
the "OBJECT" output image (CHECKIM AGE_T YPE = OB JECT) , 
that contains the detected sources separated by pixels with null val- 
ues. 10000 rectangular apertures have been randomly distributed 
over the noiseless and the original ultradeep VIMOS U-band im- 
ages. On one side, the rectangle width (W) mimic the slit width and 
on the other side, the rectangle length (L, spatial direction) repre- 
sents the inability to discriminate foreground sources closer than 
the seeing (L-2x seeing). S06 and Steidel et al. (2001) reported a 
typical seeing of > 0.8" and used a slit width of 1.2" and 1.4", 
respectively. Considering these values the simulations have been 
performed adopting rectangular sizes (WxL) of 1.2"xl.6" and 
1.4"xl.6". Since in the real case relatively bright neighbours are 
easily recognised as contaminants, we have excluded from the pos- 
sible interception sources brighter than a given U magnitude Umax 
(with Umax=24, 25 and 26). In this respect, it is important to check 
which effect has this cut in the optical magnitude domain, i.e. how 
many relatively bright sources in the optical are faint in the U band, 
and therefore recognisable as potential contaminants in the refer- 
ence band (e.g. R). The mean R magnitude of the LBGs samples 
of S06 and Steidel et al. (2001) are 23.92+0.32 and 24.33+0.38, 
respectively. Using the deep VLT/VIMOS R-band catalog of Non- 
ino et al. (2009) that reaches R-mag^29 at lcr within the 1" radius 
aperture, it turns out that with U-mag>24 (i.e. Umax-2A) the frac- 
tion of galaxies with R-mag < 24 (25) is 7%(21%). In the case of 
Umax=25 and 26 the fractions are 3%(13%) and 1%(5%) with R- 
mag < 24 (25), respectively. Therefore, in the following we refer 
to the last two cases, and in particular to the case with Umax-26 
when calculating the quantities related to the LyC contamination. 

The panels in the left- side of Figure [2] show the flux (AB) dis- 
tributions of the 10000 random positioned rectangular apertures. 
The median and the dispersion of these distributions quantify the 
contamination to the LyC estimation for the single observation, the 
dotted line shows the distribution derived from the original U-band 
image and the solid line is that of the noiseless image. The differ- 
ence between the two is clearly visible around the zero flux value, 
where in the case of the noiseless image there is a sharp peak, 
while in the real image case the peaked is around the zero value, 
but with a dispersion that is the contribution of fainter sources 
(U-mag>28.5) and background noise fluctuation (see negative side 
of the distribution). More importantly, in both cases the asymme- 
try towards brighter fluxes is the signature of pollution by inter- 
cepted sources (modulated by the Umax-26 in the case shown). 
In the following we adopt the noiseless image for the estimation 
of the LyC contamination. In order to simulate the effect in the 
average process, 10000 averages adopting the observational con- 
ditions of S06 (14 sources, W=1.2" and seeing 0.8") and Steidel 
et al. 2001 (29 sources, W=1.4" and seeing 1.0") have been per- 
formed. The right side of Figure |2] shows these distributions drawn 
from the noiseless image and adopting Umax=26, 25, 24 (solid, 
dashed and dotted lines, respectively). The medians and the one 
and two sigmas percentiles are reported in Table O Adopting the 
number counts down to 28.5, in the Steidel et al. (2001) configu- 
ration (WxL = 1.4"x2.0") and considering those realizations with 
Umax-26, the median, the lcr and 2cr lower bounds of the dis- 
tribution of the averages are 29.56, 29.03 and 28.62, respectively 
(they are [29.07, 28.48, 28.03] adopting Umax=25). Similarly, in 
the S06 conditions (WxL = 1.2"xl.6") and considering those real- 
izations with Umax-26, the median, the lcr and 2cr lower bounds 
are 30.02, 29.13 and 28.54, respectively (they are [29.62, 28.72, 



Table 2. Medians and 1, 2 cr estimations of the integrated flux due to fore- 
ground sources for the samples of S06 (14 galaxies) and Steidel et al. 2001 
(29 galaxies). The values have been derived from the noiseless images with 
number counts down to 28.5, 29.5 and 30.5. The latter has been built by 
extrapolating the observed number counts (see text for details). 



Noiseless 




N=14 






N=29 




Image 


Median 


lcr 


2(T 


Median 


lcr 


2(T 


Umax=25 


U<28.5 


29.62 


28.72 


27.92 


29.07 


28.48 


28.03 


U<29.5 


29.25 


28.54 


27.83 


28.74 


28.29 


27.90 


U<30.5 


29.13 


28.49 


27.80 


28.63 


28.21 


27.85 


Umax=26 


U<28.5 


30.02 


29.13 


28.54 


29.56 


29.03 


28.62 


U<29.5 


29.53 


28.89 


28.40 


29.08 


28.71 


28.39 


U<30.5 


29.38 


28.81 


28.35 


28.95 


28.61 


28.31 



27.92] adopting Umax=25, see Figure |2] right panels and Table |2] 
for a summary). 

It is interesting to explore how these distributions change if 
the number counts are extrapolated down to fainter flux limits 
(Umag=30.5). This has been done by fitting a straight line to the 
observed counts between U-mag 25 and 28 with a slope of 0.2 in the 
[log(N) vs. U-mag] plane (Table Q] shows the extrapolated values). 
From the real U-band image ("OBJECT" output image provided 
by SExtractor) 30 templates sources with U-mag^27 have been 
extracted. Starting from this sample the magnitude range 28.5<U- 
mag<30.5 has been populated by inserting randomly the galaxies 
appropriately dimmed according with the expected (extrapolated) 
number counts. 

The resulting medians, 1 and 2 cr of the distributions down 
to these new magnitude limits are reported in Table |2] Increasing 
the U-mag limit for the number counts, it is apparent that the me- 
dian values change more significantly than the dispersions (cr) and 
it is due to the fact that the probability to intercept (faint) sources 
significantly increases (the probability to intercept a contaminat- 
ing source is >20%). Conversely, the variation of the scatter (cr) is 
lower, since it is mainly sensitive to the brighter contaminants, as 
expected when an average is computed. It is also worth noting that 
passing from U-mag 28.5, 29.5 and 30.5, the relative variations of 
the medians and the dispersions decrease, suggesting that a further 
extrapolation (e.g. down to U-mag 31.5) do not change significantly 
the present results. 

It is now possible to compare the observed resid- 
ual flux reported in literature with the present simula- 
tions. Steidel et al. (2001) measured a residual flux in 
their composite spectrum of 29 LBGs of m AB -21 A (AB) 
(/900(AB)=4.02xlO- 31 erg sec^cm^Hz' 1 ) and no direct 
detection was identified. They reported an observed ratio 
(/1500//900)obs of 17.7+3.8 that corresponds to a contrast of 3.1 
magnitudes in the AB scale. S06 using a slit width of 1.2" and the 
average over the 14 sources measured 28.33 (AB) also shown in 
the Figure |2] (adopting their observed ratio /1500//900=58 and 
the average magnitude of the sample i?=23.92). Considering the 
simulation with the noiseless image and number counts down to 
U-mag=28.5, it turns out that the LyC estimations by Steidel et 
al. (2001) and S06 are significant at >2cr level from the expected 
median foreground contamination. Adopting the resulting median 
flux values derived from the simulations it is now possible to 
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correct the LyC value observed by S06. Assuming the case of 
the noiseless image with number counts down to 30.5 and Umax 
= 25 and 26, it turns out that the corrected residual /900 flux 
would be ~ 1.9 (29.04 AB) and 1.6 (28.85 AB) times smaller 
than the value observed (28.33 AB), respectively. In the case of 
Steidel et al. (2001) the average residual /900 flux (magnitude 
27.4) is brighter than S06, so the correction is less important. The 
corrected residual /900 flux would be ~ 1.5 (27.82 AB) and 1.3 
(27.70 AB) times smaller, adopting the medians calculated from 
the noiseless images down to U-mag 30.5 with Umax = 25 and 26, 
respectively. This would also decrease the average relative escape 
fraction (< f esc ,rei >) by the same factors. 

A way to further investigate the contamination occurrence due 
to foreground blue sources on z ~ 3-4 galaxies is to use high- 
resolution, deep and multicolour images as those obtained from the 
GOODS project. This is the argument of the next section. 



3 U-BAND DETECTED LEGS AT Z>3.48 IN THE 
GOODS-S 

An observational test can be performed directly on the GOODS - 
S area by cross-correlating the VLT/VIMOS U-band catalog with 
the list of available B-band dropout galaxies observed during 
the VLT/FORS2 campaign with known spectroscopic redshift 
(Vanzella et al. (2008, 2009a)) and for which the Lyman limit is red- 
der than the filter transmission. The astrometric solution between 
the U and the ACS catalogs gives an r.m.s. in RA and DEC lower 
than 0.1" on average (Nonino et al. 2009). Adopting a radius of 
0.5", 1.0" and 1.3" in the matching, it turns out that 1, 4 and 11 
galaxies out of 36 with redshift between 3.47-4.3 have an U-band 
detection with U-mag in the range 26.1-29.40 the maximum of the 
probability distribution peaks at 1 and 4 galaxies contaminated out 
of 36. Therefore the numbers we obtain matching the U-band and 
the ACS catalog are fully consistent with the random occurrence of 
foreground superpositions (see Table [3]). 

In the following we consider the matching within 1" (since 
typically ground-based imaging is performed with seeing better 
than 1"). In this case the four LBGs that lie within 1.0" are: 
(1) GDS J033225. 16-274852.6 (U=26.57), (2) GDS J033220.97- 
275022.3 (U=27.17), (3) GDS J033226.28-275245.7 (U=28.47) 
and (4) GDS J033236.83-274558.0 (U=28.63) (see Figure[3]). The 
magnitudes of the involved U-band detected sources range within 
26.57<U<28.63, corresponding to a percentage of the observed 
continuum flux of the LBGs in the z 7 75 band (~ 1600- 1700 A rest- 
frame) of 2-10%, a quantity that resemble the observed LyC val- 
ues often reported in literature. Apart from the brighter case (1) in 
which the contamination comes from the relatively bright galaxy 
close to the LBG (GDS J033225.09-274852.6, z 850 = 22.61, B 435 = 
26.0), in the other three cases the depth of the GOODS (and HUDF) 



2 See Steidel et al. (2001) or S06 Eq. 1 for the definition. 

3 At redshift larger than 3.4 the 912A limit is redder than the VIMOS U- 
band response (Nonino et al. 2009). It is also known that one of the four 
quadrants of the VIMOS U-band image shows a read leak at i~4850A, 
however dedicated simulations show that it has a negligible effect (^0.1%) 
(see Nonino et al. (2009), appendix C). Moreover the quadrant is not affect- 
ing the cases we are discussing here, especially that on the HUDF. In the 
following we consider the cases with U-mag < 28.5, for which the detection 
is more reliable (>3<x). Given the sample of LBGs considered (N=36) and 
adopting the number counts down to U<28.5, 
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Figure 3. BVi composite ACS colour images (top), the VLT/VIMOS U 
band images (bottom) and the VLT/VIMOS R band images (bottom) of the 
four cases (out of 36 LBGs) that satisfy the matching between the U-band 
and the ACS v2.0 catalogs, adopting a radius of 1.0" (the cutouts of 10" size 
show several blue compact sources in the blank field). From left to right the 
U magnitude increases and the arrows indicate the position of the LBGs. 
The upper-right case is shown with more detail in Figure |4] (dotted square). 



Table 3. Probabilities % P(K) and ¥(>K)) that "K" and "at least K" galax- 
ies out of 36 (AO in the GOODS-S field used to check the occurrence of 
foreground contamination within 0.5" and 1.0" are reported (assuming U- 
band number counts down magnitude 28.5). The number of cases observed 
(1/36 with 0.5" and 4/36 with 1.0" radius) are highlighted in bold face. 







?(>K) 


P(K) 


P(K) 


K 


0.5" radius 


1.0" radius 


0.5" radius 


1.0" radius 





100.0 


100.0 


31.6 


0.8 


1 


68.4 


99.2 


37.0 


4.1 


2 


31.4 


95.2 


21.1 


10.2 


3 


10.4 


84.9 


7.8 


16.8 


4 


2.6 


68.2 


2.1 


19.9 


5 


0.5 


48.2 


0.4 


18.4 



area allow to visually appreciate the presence of at least one com- 
pact blue source (not distinguishable in the ground-based R band) 
shifted to respect the targeted LBG and within 1" distance. In par- 
ticular in the case (4), the fainter in the U-band among the four 
discussed, the spatial offset is < 0.5", and sits in the HUDF. In the 
next section we report on this (extreme) system as an example of 
LyC contamination. 

3.1 An exemplary faint foreground contamination at z=3.797 
in the HUDF 

An example of false detection of escape radiation of one of the B- 
band dropout galaxies discussed above is shown in the right panel 
of Figure |4] The clear detection in the U band is shown (U mag- 
nitude 28.63+0.2) at the position of the galaxy GDS J033236.83- 
274558.0 at redshift 3.797. In the left panel of the same figure the 
colour image extracted from the HUDF is shown. It is clearly vis- 
ible a blue compact source superimposed on the z=3.797 galaxy 
slightly offset by <0.4" (this blue source is visible also in the 
ACS/B435 and V606 bands). We note that the LBG and the blue 
source have been detected as a single galaxy in the HUDF catalog 
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Blue source 






ACS HUDF zoomed 



ACS HUDF VIMOS U-band 



(c) 



Figure 4. Left: Detail of ACS HUDF colour image of the B-band dropout 
galaxy (case 4 in Figure at redshift 3.797. A blue compact source at 
-0.4" from the center of the LBG is indicated, superimposed to the high 
redshift galaxy. Middle: 10" side field around the LBG. The yellow dotted 
squared is shown in the left panel. Other three compact U-detected sources 
are marked, a), b) and c). Right: The VLT/VIMOS U-band image of the 
middle panel (Nonino et al. 2009). The detection at the position of the high 
redshift galaxy is clear (5<x, solid circle), in which the presence of a non- 
resolved (in the U-band) blue compact source is evident in the HUDF image 
(left side). 



by Coe et al. (2006). If the blue spot is not a contaminant, then it 
would be an ultradeep, high resolution "morphologically detected" 
LyC case at high redshift. It is not clear what would be the morpho- 
logical appearance of galaxies showing LyC escape photons at A < 
912A, however the system can be easily explained by the presence 
of a foreground source rather than a shifted "hole" where Lyman 
photons can escape. Therefore we favour the simpler interpretation 
of a superposition case. It is worth noting that this system would 
appear as a single source if observed with lower spatial resolution, 
as in the R band or deep U image (seeing 0.8", where presumably 
the blue source is entirely contributing to the U-band detection). A 
deep spectroscopic or narrow band observation of this system with 
a seeing larger than 0.5" would detect a signal below 912A erro- 
neously ascribed to the source at higher redshift. In particular, given 
the U mag of 28.63 and 1775 = 24.67 the polluting source mimic an 
observed (f 1500/ f 900) bs ~ 38. This value is at least three times 
larger than the current LyC detections (S06 and 109), and compa- 
rable to the reported estimates of stacked spectra, i.e. 17.7 and 58 
reported by Steidel et al. (2001) and S06, respectively. 



4 CONCLUDING REMARKS 

Depending on the seeing conditions and/or the aperture sizes 
adopted to carry out the photometry/spectroscopy the contamina- 
tion probability by foreground blue sources on the estimation of 
the ionising radiation may not be negligible. In particular in the 
cases reported by S06 and 109 it turns out that the probability of 
foreground contamination for at least half of the direct LyC detec- 
tions is ~ 30% and 50%, respectively. Among the direct LyC detec- 
tions in 109 there are cases of spatially offset emission (in the case 
of LBGs) and cases of extremely high escape fraction (even larger 
than unity). We argue that, other than a possible complex physical 
explanations (e.g. Inoue 2009), these suspicious systems should be 
better investigated with higher resolution and multi-band imaging, 
since they are statistically compatible with superposition effects. 

Also in the case of statistical LyC detection the foreground 
contamination plays a role in the global result. A stacking process 
increases the signal to noise of the possible true LyC emission, but 
at the same time a bias by foreground faint sources (not detected in- 
dividually) may arise. The issue has been addressed by computing 
Monte-Carlo simulations from which, taking advantage of the ul- 
tradeep VLT/VIMOS U-band number counts (Nonino et al. 2009), 



we calculated the median contribution of foreground sources to the 
residual flux measure below the Lyman limit. From these calcula- 
tions it turns out that the current estimations of LyC emission from 
stacked spectra at redshift ~3 are beyond 2<x from the expected 
median foreground contribution. 

The simple foreground contamination may not entirely explain 
the cases of direct escape fraction observed or the limits derived 
from the stacking, nevertheless the calculations and the examples 
here reported demonstrate the need for caution in LyC measure- 
ments at high redshift. This systematic error has to be considered 
when attempting to measure the LyC escape fraction and its evolu- 
tion with redshift and the contribution of galaxies to the UV back- 
ground. For example, from simulations it turns out that the cor- 
rected observed /900 flux reported by S06 would be ~ 1.6-1.9 
times smaller than the reported value. Differently, the correction for 
the Steidel et al. (2001) estimation is less effective (a factor ~ 1.3), 
since their residual flux is ~ 1 magnitude brighter than S06. Conse- 
quently the relative escape fractions are dimmed by the same fac- 
tors. It is worth noting that in stacked images of tens and hundreds 
galaxies at redshift ~ 1 there is currently no direct detection (Cowie 
et al. (2009), Siana et al. (2007, 2010), Malkan et al. (2003)). On 
the one hand, it may reflect indirectly the redshift distribution of 
the (faint U-mag>26) blue compact sources that populate the ultra- 
deep U-band image. Since not evident contaminations have been 
currently measured at redshift ~ 1, the bulk of the (faint) blue com- 
pact sources would be at redshift > 1. On the other hand, being 
these objects U-band detected, their typical redshift is most prob- 
ably smaller than 3 (i.e. they are not U-band dropouts). For the 
brighter part of the distribution (U-mag<26.5), a further constraint 
on the redshift comes from the MUSIC photometric redshift cata- 
log (Grazian et al. 2006), from which they span the interval 0<z<3. 

Conversely, in the case of LyC observations at z>3 the con- 
tribution of such blue galaxies as contaminants largely increases. 
In general, the fraction of galaxies contaminated ranges between 
~3 and ~15%, varying from seeing 0.5" to 1.0", respectively, and 
assuming number counts down to U-mag~28.5. On the one hand 
simulations predict an evolution (increase) of the escape fraction 
with increasing redshift (e.g. Razoumov & Sommer-Larsen 2007), 
on the other hand also the contamination by foreground sources it 
is expected to increase, with a more significant effect for fainter 
LBGs (R > 24). Therefore, if the contamination is not taken into 
account properly, the observational constraints on the evolution of 
the escape fraction may be biased, particularly at faint limits where 
the population of blue compact sources increases in number den- 
sity. In order to produce a robust LyC measurements at z>3 high 
spatial resolution, multi-band and ultradeep imaging are necessary 
to exclude the spurious cases. 
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